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The variational method allows one to study the mixing of interpolators with different chiral trans- 
formation properties in the nonperturbatively determined physical state. It is then possible to 
define and calculate in a gauge-invariant manner the chiral as well as the partial wave content 
of the quark-antiquark component of a meson in the infrared, where mass is generated. Using a 
unitary transformation from the chiral basis to the 2s+l Lj basis one may extract the partial wave 
content of a meson. We present results for the p- and p '-mesons using a simulation with Nf = 2 
dynamical quarks, all for lattice spacings close to 0.15 fm. Our results indicate a strong chiral 
symmetry breaking in the p state and its simple 3 5i-wave composition in the infrared. For the p'- 
meson we find a small chiral symmetry breaking in the infrared as well as a leading contribution 
of the 3 Di partial wave, which is contradictory to the quark model. 
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1. Introduction 

A central question in QCD is the mass generation mechanism and its interconnection with con- 
finement and chiral symmetry breaking. The angular momentum generation of hadrons is another 
related question. Chiral symmetry is dynamically broken in the QCD vacuum which is evidenced 
by the absence of parity doublets in the low-lying hadron spectrum and by the existence of the 
pion as a pseudo Goldstone boson. It follows from the trace anomaly that the hadron mass al- 
most completely consists of the energy of the quantized gluonic field. However, this statement 
tells us nothing about mechanism of mass generation. It is widely believed that the spontaneous 
breaking of chiral symmetry, i.e., the quark condensate of the vacuum, is prominently responsible 
for the mass generation of hadrons like the nucleon or p -meson. This was observed in various 
microscopical models and in the QCD sum rule approach [jl|, ||]. 

Chiral symmetry breaking in the vacuum explains a phenomenological success of the quark 
model, at least for the ground states. Namely, almost massless light quarks acquire their effective 
(dynamical, constituent) mass at low momenta via their coupling with the quark condensate. This 
large mass renders the problem effectively non-relativistic and the ground state p is a 3 S\ state in 
the quark model language [^, H]. Traditionally the excitation spectrum is also described within the 
quark model and the first excited state of the p-meson, p(1450), is believed to be the first radial 
excitation, i.e., a 3 S\ state 1^,^]. 

At the same time there are indications that chiral symmetry is effectively restored in the high- 
lying spectrum [Q, ^]. This would imply that the mass generation mechanism in highly excited 
hadrons is different and the quark condensate of the vacuum is of little importance. It would also 
imply that the constituent quark model language is inadequate for excited hadrons. To resolve the 
issue one needs direct information about the hadron structure, which can be supplied in ab initio 
lattice simulations. Here we present a way to reconstruct in dynamical simulations a chiral as well 
as an angular momentum decomposition of the leading quark-antiquark component of mesons at 
physical, infrared scale. 

The variational method |{7|] provides a tool to study the hadron wave function in lattice QCD 
calculations. One uses a set of interpolators {0\ ,02,--, On}, which have the quantum numbers 
of the state of interest, and computes the cross-correlation matrix, 



One solves a generalized eigenvalue problem; assuming that the set of interpolators {£?,} is com- 
plete enough, the wave function is related to the eigenvectors obtained. We are interested in the 
reconstruction of the leading quark-antiquark component of the low lying mesons. Therefore we 
need interpolators that allow us to define such a component in a unique way. 

In [||] a classification of all non-exotic quark-antiquark states (interpolators) in the light meson 
sector according to the transformation properties with respect to the SU (2)l X SU (2)r and U(1)a 
was presented. If no explicit excitation of the gluonic field with non-vacuum quantum numbers is 
present, this basis is a complete one for a quark-antiquark system and we can define and investi- 
gate chiral symmetry breaking. Namely, one can reconstruct a decomposition for a given meson 
in terms of different representations of the chiral group by diagonalizing the cross-correlation ma- 
trix from < \l.l\ ). The eigenvectors describe the quark-antiquark content in terms of different chiral 
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representations. If we observe components with different transformation properties in terms of 
SU (2)l x SU (2)r and U(1)a, then we conclude that chiral symmetry is broken in that state. 

In order to establish a connection to the quark model, it is interesting to reconstruct the meson 
composition in terms of the 2S+i Lj basis, where J = L + S are the standard angular momenta in the 
two-body system. Such a decomposition of the leading quark-antiquark component in terms of the 
2S+l Lj basis in the infrared, i.e., where the hadron mass is generated, would tell us to which degree 
the quark model language is adequate for a given state. 

The 2S+l Lj angular momentum basis and the chiral basis are both complete for a two-particle 
system. They are connected to each other via a unitary transformation [§]. Each state of the 
chiral basis can be uniquely represented in terms of the 2S+l Lj states. Then, diagonalizing the 
cross-correlation matrix, built from interpolators with definite chiral transformation properties, one 
can obtain the partial wave decomposition of the leading Fock component. This method can in 
principle be applied to any meson, here we use as an example the vector meson p and its first 
excitation p (1450) [§, |l(J. 

2. Chiral classification and the transformation to the angular momentum basis 

The classification of the quark-antiquark states and interpolators with respect to representations 
of SU(2)l x SU(2)r was done in [Q]. We are interested in the quark-antiquark component of 
the ground state p -meson and its first excitation. There are two possible chiral representations 
that are compatible with the quantum numbers of the p -meson, which have drastically different 
chiral transformation properties. Assuming that chiral symmetry is not broken, then one has two 
independent states. The first state is |(0, 1) © (1,0); 11 ); it can be created from the vacuum by 
the standard vector current 

O v = qfxq. (2.1) 

Its chiral partner is the a\ meson. The other state is |(1 /2, 1 /2)&; 1 1 ), which can be created by 
the pseudotensor operator, 

O T =qa 0i rq, (2.2) 

and its chiral partner is the h\ meson. Here, f denotes the vector of isospin Pauli matrices. 

Chiral symmetry breaking in the state implies that the state should be a mixture of both repre- 
sentations. If it is a superposition of both representations with approximately equal weights, then 
the chiral symmetry is maximally violated in this state. If, on the contrary, one of the represen- 
tations strongly dominates over the other representation, one could speak about effective chiral 
restoration in this state. 

These chiral representations can be transferred into the 2S+l Lj basis, using the unitary trans- 
formation [||] 



where U is given by 
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Thus, using the interpolators Oy and Or from (2.1) and (2.2) for the diagonalization of the 
cross-correlation matrix, we are able to reconstruct the partial wave content of the leading qq-Fock 
component of the p -meson. 



3. Reconstruction of the wave function using the variational method 

We briefly want to discuss the basic features of the variational method [^] and how to analyse the 
decomposition of the p -mesons. The time propagation properties of the normalized physical states 
\n) are given by 

(«(0K0)) = <5„ m e- £ " ; . (3.1) 

The lattice interpolators 0, are typically not normalized and are projected to zero spatial momen- 
tum. The cross-correlation matrix from Eq. ( LI ) can be written as 

Qj(t) = (O i (t)O}(0)) = ^af*^ , (3.2) 

n 

where the coefficients give us the overlap of the physical state \n) with the lattice interpolator 

Ou 

af ] = (0|O,-|b) . (3.3) 

The two chiral representations (0, 1) (1,0) and (1/2, 1 /2)t, form a complete and orthogonal basis 
(with respect to the chiral group) for p -mesons. Consequently, using the variational method we are 
able to study the mixing of the two representations in both p and p' states. 

Following the lines of [|J] one can show that the ratio of couplings can be obtained as 

a<r\t)_ Qj(t)u?(t) 

4 n) (o c kj {t)uf{ t ) ' 

The ratio of the vector to pseudotensor couplings, ay /a^ , tells us about the chiral symmetry- 
breaking in the states n = p,p'. 



4. Defining the resolution scale 

If we probe the hadron structure with the local interpolators, then we study the hadron decompo- 
sition at the scale fixed by the lattice spacing a. For a reasonably small a this scale is close to the 
ultraviolet scale. However, we are interested in the hadron content at the infrared scales, where 
mass is generated. For this purpose we cannot use a large a, because matching with the continuum 
QCD will be lost. Given a fixed, reasonably small lattice spacing a a large infrared scale R (i.e., 
small resolution scale l/R) can be achieved by gauge-invariant smearing of the point-like interpo- 
lators. We smear the quark field (sources) in spatial directions over the size R in physical units, 
such that R/a » 1 . Then, even in the continuum limit a — > we probe the hadron content at the 
infrared scale fixed by R. Such a definition of the resolution is similar to the experimental one, 
where an external probe is sensitive only to quark fields (it is blind to gluonic fields) at a resolution 
that is determined by the momentum transfer in spatial directions. 
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Set 


Plw 




#conf 


a [fm] 


m n [MeV] 


m p [MeV] 


A 


4.70 


-0.050 


200 


0.1507(17) 


526(7) 


911(11) 


B 


4.65 


-0.060 


300 


0.1500(12) 


469(5) 


870(10) 


C 


4.58 


-0.077 


300 


0.1440(12) 


323(5) 


795(15) 



Table 1: Details of the lattice simulation: the leading value fi LW of the gauge coupling, the bare mass pa- 
rameter a mo of the CI action, the number of analyzed configurations, the lattice spacing a, the pseudoscalar 
mass m K , and the vector meson mass m p (cf. Jl3[ ] for more details). 



The smearing procedure itself is done using Jacobi smearing Qll|]. A smearing operator M acts 
on a point-like source So, 

N 

S = MS , M=Y 4 { kH Y- (4- 1 ) 



n=0 



The hopping term H is given by 



3 

rt 



H= £ [ Uk (x,t)5 x+h + Ul(x-k,t)S x _ lY ] . (4.2) 



',y 
k=\ 



It creates approximately a Gaussian profile of the width R for each quark field of the smeared 
interpolator. 

5. Lattice simulation details and results 



Like in our previous analysis of excited hadrons [112J, 13], we use the Liischer-Weisz gauge action 



| ]14p and the Chirally Improved (CI) Dirac operator, which has better chiral properties than the 
Wilson Dirac operator [O]. For this study three sets of dynamical gauge configurations, all for 
lattice size of 16 3 x 32, including two mass-degenerate light sea quarks are used (for details see 
Tab.!). 

We include in our cross-correlation matrix the four interpolators 

0\ = u n fd n , 2 = u w fd w , (5. 1) 

O3 = u n i f d n , 4 = u w i Y d w . (5.2) 

Here f denotes one of the spatial Dirac matrices and / the y-matrix in (Euclidean) time direction. 
The subscripts n and w (for narrow and wide) denote the two smearing widths, R w 0.34 fm and 
0.67 fm, respectively. With these interpolators we are able to extract both the ground state mass 
and the mass of the first excited state of the p-meson, see the l.h.s. of Fig. [jj 

On the r.h.s. of Fig. [T] we show the /?-dependence of the ratio from Eq. (^4) for the case 



ay jaj both for the ground state p-meson and its first excited state. This ratio of the vector to the 
pseudotensor coupling to the states shows us their decomposition in terms of the (0, 1 ) © ( 1 , 0) 
and (1/2, l/2)f, representations. For the ground state at the largest value of R 0.67 fm this 
ratio is approximately 1.2, i.e., we see a strong mixture of the two representations in the wave 



function of the ground state p-meson. Inverting the unitarity transformation from Eq. (2.3) results 
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Figure 1: L.h.s.: The vector meson mass my is plotted against m| for all three sets. Black circles represent 
the ground state, p, and red squares represent the first excitation, p'. The experimental values are depicted as 
magenta crosses with decay width indicated. R.h.s.: The ratio ay far is plotted against the smearing width 
R for all three data sets. Black circles represent the ground state and red squares the first excitation. Broken 
lines are drawn only to guide the eye (color online). 



in the fact that the vector meson is predominantly a 3 5i state with a tiny admixture of a 3 D\ wave, 
0.997 | 3 Si) -0.073| 3 Di). This result indicates that the p(770) at the scale fixed by the meson size 
is approximately a 3 Si state - in agreement with the quark model language. 

However, the situation changes for the first excited state, p' = p(1450). In this case a strong 
dependence of the ratio on the infrared scale is observed. Extrapolating the results to the scale of the 
p' size, R ~ 0.8 — 1 fm, one expects a significant contribution from the (1/2, 1/2)/, representation 
and a contribution of the other representation is suppressed. This indicates a smooth onset of 
effective chiral restoration. 

The interpretation is as follows. From the conformal symmetry of QCD one expects that in 
the deep ultraviolet the pseudotensor interpolator decouples from the p -mesons. This can also 
be seen from the non-vanishing anomalous dimension of the pseudotensor operator, implying its 
decoupling in the ultraviolet limit. Thus, the ratio ay / aj must increase for small R. At large R 
the ratio determines a degree of chiral symmetry breaking in the infrared region, where mass is 
generated. 

In the p (770) meson chiral symmetry is strongly broken since this state is a strong mixture 
of (0,1) ©(1,0) and (1/2,1/2)/, with approximately equal weights. Consequently, its "would- 
be chiral partners" have a much larger mass: «i(1260) and /zi ( 1 170). To these low lying states 
we cannot assign any chiral representation. For the p(1450) the contribution from (1/2, 1/2)/, is 
much bigger than the contribution of the other representation. One then predicts that in the same 
energy region there must exist an h\ (and not an a{) meson as an approximate chiral partner of 
p(1450). And in fact there is a state h\ (1380) and no a\ state in the same energy region. The second 
excited p -meson, the p(1700), should then be dominated by the representation (0, 1) © (1,0). This 
assumption is favored by the existence of the «i(1640) state. There is no room for this ai(1640) 
meson within the the quark model @|. 

Although we do not have the precise value of the ratio ay/aj for p(1450) at large R ~ 0.8 — 
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lfm, it is indicative that this value is small. Then we are able to give a qualitative estimate for the 
angular momentum content of the p ( 1450) in the infrared. Assuming a vanishing ratio the state 
would have the following partial wave content, 

] jlf Sl )-^\ 3 D l ). (5.3) 

This shows a leading contribution of the 3 Di wave. Possible small deviations of the ratio from zero 
do not change this qualitative conclusion. This result is inconsistent with p' to be a radial excitation 
of the ground state p-meson, i.e., an 3 Si state, as predicted by the quark model 

Acknowledgments 

We thank G. Engel and C. Gattringer for discussions. L.Ya.G. and M.L. acknowledge support 
of the "Fonds zur Forderung der Wissenschaftlichen Forschung" (P21970-N16, DK W1203-N08) 
and DFG project SFB/TR-55, respectively. The calculations have been performed on the SGI Altix 
4700 of the Leibniz-Rechenzentrum Munich and on local clusters at the ZID at the University of 
Graz. 

References 

[1] M. A. Shifman, A. I. Vainshtein, V. I. Zakharov, Nucl. Phys. B 147, 385 (1979). 
[2] B. L. Ioffe, Nucl. Phys. B 188, 317 (1981); erratum 191, 591 (1981). 

[3] C. Amsler, T. DeGrand, B. Krusche, "Quark Model", p 172 ff in Particle Data Group, Phys. Lett. B 
667, 1 (2008). 

[4] S. Godfrey, N. Isgur, Phys. Rev. D 32, 189 (1985). 

[5] L. Y. Glozman, Phys. Lett. B 587, 69 (2004); Phys. Rep. 444, 1 (2007). 

[6] L. Y. Glozman, Phys. Rev. Lett. 99, 191602 (2007). 

[7] C. Michael, Nucl. Phys. B 259, 58 (1985); M. Luscher, U. Wolff, Nucl. Phys. B 339, 222 (1990). 

[8] L. Y. Glozman, A. V. Nefediev, Phys. Rev. D 76, 096004 (2007); ibid. 80, 057901 (2009). 

[9] L. Y. Glozman, C. B. Lang, M. Limmer, Phys. Rev. Lett. 103, 121601 (2009); Few Body Syst. 47, 91 
(2010). 

[10] L. Y. Glozman, C. B. Lang, M. Limmer, arXiv: 1007. 1346 [hepdat], Phys. Rev. D, in print. 

[11] S. Giisken et al, Phys. Lett. B 227, 266 (1989); C. Best et al., Phys. Rev. D 56, 2743 (1997). 

[12] T. Burch et al., Phys. Rev. D 70, 054502 (2004); T. Burch et al., Phys. Rev. D 73, 094505 (2006); 
C. Gattringer et al., Phys. Rev. D 78, 034501 (2008). 

[13] C. Gattringer et al., Phys. Rev. D 79, 054501 (2009); G. P. Engel et al., Phys. Rev. D 82, 034505 
(2010). 

[14] M. Luscher, P. Weisz, Commun. Math. Phys. 97, 59 (1985). 

[15] C. Gattringer, Phys. Rev. D 63, 1 14501 (2001); C. Gattringer, I. Hip, C. B. Lang, Nucl. Phys. B 597, 
451 (2001). 



7 



